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I  INTRODUCTION 

For  a  number  of  years  the  personnel  at  AeroChem  have  been  protagonists  of 
the  ionic  mechanism  of  soot  formation  in  flames.  In  this  theory,  it  is  as¬ 
sumed  that  chemi-ions  grow  very  rapidly  through  a  series  of  ion-molecule  reac¬ 
tions  to  produce  large  ions  which  are  neutralized  to  produce  either  very  small 
incipient  soot  particles,  or  very  large  neutral  molecules  which  rapidly  grow 
to  form  the  incipient  soot  particles.  We  have,  in  fact,  previously  demon¬ 
strated  that  there  is  no  sharp  demarcation  between  large  molecules  and  small 
soot  particles. 1 After  much  study  by  many  people,  an  understanding  of  the 
soot  nucleation  process  still  represents  one  of  the  more  significant  problems 
challenging  the  combustion  community.  A  large  part  of  the  reason  for  this  is 
that  the  necessary  data  required  to  resolve  the  many  questions  about  the  ini¬ 
tial  stages  of  soot  formation — where  the  major  questions  apply — are  very  dif¬ 
ficult  to  obtain  experimentally;  global  measurements  of  soot  particles  do  not 
give  the  necessary  information.  The  objective  of  our  program  has  been  to  make 
the  necessary  detailed  ion  (both  total  and  individual)  concentration  profile 
measurements  in  the  same  flames  in  which  Howard  and  associates  at  MIT3**  and 
others,  e. g. ,  Delfau  and  associates7*9  in  France  and  Homann  and  associates  in 
Germany13**3  are  making  complementary  measurements.  Our  measurements  include: 
mass  spectrometer  measurements  for  individual  ion  concentrations;  Langmuir 
probe  measurements  for  total  ion  concentrat ions;  and  temperature  compensated 
thermocouple  measurements  for  temperature  profiles.  To  assure  that  the 
flames  we  study  are  exactly  comparable  to  those  used  by  Howard  and  associates 
at  MIT,  we  have  duplicated  their  burner. 

During  this  reporting  period  we  were  honored  by  an  invitation  to  present 
our  work  on  the  ionic  mechanism  of  soot  formation  at  the  Gordon  Conference  on 
Fuels.  As  a  theme  we  used:  "Why  I  Believe  In  The  Ionic  Mechanism  Of  Soot 
Formation  In  Flames”.  The  evidence  was  reviewed  there  in  detail;  for  the 
record  it  is  summarized  in  Table  I. 

The  work  statement  for  this  program  contained  two  phases.  Phase  I  expei — 
iments  are  to  be  done  in  well  studied  flames  of  benzene-oxygen  and  acetylene- 
oxygen  and  they  include  generation  of  temperature  profiles,  Langmuir  probe 
curves  of  total  ion  concentrat ion,  and  mass  spectrometer  profiles  of  individ¬ 
ual  ions.  Phase  II  involves  model  development  of  the  individual  processes 
involved  in  incipient  soot  formation,  e. g.,  rates  of  ion  and  particle  forma¬ 
tion,  rates  of  production  of  large  molecular  ions,  the  temperature  of  the 
growing  particle  and  its  effect  on  thermal  ionization,  thermal  ionization  of 
large  molecules,  and  ion  molecule  rate  coefficients.  The  quantitative  models 
developed  in  this  program  for  each  of  the  steps  will  form  the  basis  for  a 
detailed  overall  quantitative  model  of  the  total  process  of  soot  formation. 
Considerably  more  time  was  consumed  in  the  calibration  of  the  mass  scale  than 
anticipated  but  this  was  necessary  in  order  to  treat  the  problem  quantita¬ 
tively.  This  was  done  at  the  expense  of  other  experimental  tasks. 


In  this  report  we  summarize  the  results  on  this  program  during  this  report 
period.  Result  which  have  been  published  will  be  only  briefly  summarized; 
other  results  will  be  covered  in  more  detail. 


II.  RESEARCH  STATUS 

A.  APPARATUS 

The  apparatus  used  in  the  mass  spectrometric  studies  reported  here  con¬ 
sists  of  a  multitubular  burner  supported  on  a  vertical  translation  stage  which 
feeds  through  a  vacuum  seal  into  a  cylindrical  flame  chamber  pumped  by  a  large 
(140  L's**)  mechanical  pump.  The  top  of  the  chamber  is  fitted  with  a  water — 
cooled  plate  supporting  a  metal  sampling  cone  which  permits  the  sampling  of 
flame  gases  through  an  (0.013-0.025  cm)  orifice  into  the  first  of  two  differ — 
entially  pumped  vacuum  chambers.  This  chamber,  maintained  at  less  than  0.01 
Pa  pressure,  contains  a  series  of  electrostatic  lenses  which  focus  and  direct 
the  flame  ion  beam  into  the  second  vacuum  chamber  containing  a  quadrupole  mass 
filter,  maintained  at  pressures  below  0.001  Pa  by  a  15  cm  (6  in.)  diffusion 
pump.  The  quadrupole  used  in  this  work  was  constructed  at  AeroChem  about  20 
years  ago.  In  this  study  it  was  powered  by  a  350  kHz  rf  supply  providing  a 
useful  mass  range  up  to  nearly  600  amu.  Ions  which  pass  through  the  mass 
filter  strike  the  cathode  of  an  electron  multiplier  and  the  amplified  current 
is  recorded  on  an  XY  recorder  as  a  function  of  mass  (fixed  burner  position)  or 
as  a  function  of  burner  position  (fixed  ion  mass). 

The  copper  burner  used  by  Bittner  and  Howard*  at  MIT  was  duplicated*4  so 
we  could  compare  our  ion  profiles  with  their  profiles  of  neutral  species  in 
the  same  flames. 

During  this  report  period  we  acquired  a  digitizer  pad  so  that  we  could 
transfer  mass  spectrometer  data  from  recorder  charts  to  a  computer  for  proces¬ 
sing. 


B.  MASS  SPECTROMETER  MASS  CALIBRATION 

Because  we  extended  the  mass  spectrometer  range  to  fairly  high  mass 
ranges,  approaching  600  amu,  it  was  necessary  to  calibrate  the  mass  range. 

This  calibration  was,  in  fact,  complicated  by  the  high  range  because  no  good 
standard  ions  are  available  in  this  range.  Normally,  the  introduction  into 
the  flame  of  easily  ionized  alkali  metal  salts  provides  unambiguous  mass  mark¬ 
ers  in  the  hot  flame  gases.  In  this  case,  the  salts  of  potassium  (39  and  41 
amu),  rubidium  (65  and  87  amu),  cesium  (133  amu),  and  lead  (206  and  208  amu) 
were  useful  only  in  the  low  mass  region.  Due  to  strong  nonlinearity  of  the 
mass  scale  below  50  amu  and  more  subtle  nonlinearity  in  other  ranges,  reliable 
extrapolations  to  600  amu  were  not  feasible  even  with  the  use  of  the  above 
salts  and  the  many  known  hydrocarbon  masses  below  about  250  amu  previously 
observed  in  rich  acetylene-oxygen  flames  with  greater  mass  resolution.7’**  In 
order  to  identify  the  masses  of  the  many  heavy  ions  observed  beyond  200  amu, 


we  recorded  and  compared  the  spectra  from  fuel  rich  deuterated  acetylene 
flames.  In  order  to  conserve  the  isotopic  fuel,  a  smaller  burner  was  selected 
for  the  calibration  experiments.  This  burner  was  similar  to  the  stainless 
steel  burner  used  in  previous  studies1** 1 7  except  that  the  central  tubes  were 
fed  gases  separately  from  the  outer  tubes  which  gave  a  2. 64  cm  central  burner 
surrounded  by  a  4.70  cm  o. d.  annular  burner.  The  annular  flame  (C2H2/02) 
provided  a  simple  and  rapid  way  to  switch  gases  in  the  central  burner  and 
quickly  reignite  the  test  flames  (C2H2/02  or  C2D2/02).  Once  the  central  flame 
became  stable,  the  annular  flame  was  extinquished  to  avoid  diffusional  mixing 
of  protonated  and  deuterated  species.  Attempts  to  stabilize  SO  cm/s  (unburned 
gas  velocity)  flames  at  2.67  kPa  were  fruitless  on  this  small  diameter  burner, 
so  a  velocity  of  about  60  cm/s  was  selected  which  provided  stable  flames  at  an 
equivalence  ratio  of  2.5  to  2.6.  These  flames  are  rich  in  ions  throughout  the 
mass  range  of  the  spectrometer. 

The  C2D2,  was  produced  from  the  reaction  of  D20  (99.6%  D,  Stohler  Iso¬ 
tope  Chemicals)  with  calcium  carbide.  The  C2D2  passed  through  a  drying  tube 
into  a  5  L  storage  vessel.  After  this  vessel  was  filled  to  1  atm,  the  genera¬ 
tion  was  terminated  and  the  vessel  connected  to  a  calibrated  sonic  orifice 
flowmeter.  The  C2D2  was  fed  to  the  burner  at  a  constant  rate  by  metering 
nitrogen  into  a  balloon  in  the  storage  vessel  to  maintain  a  constant  pressure 
(about  90  kPa  absolute)  upstream  of  the  0.025  cm  orifice.  The  downstream 
pressure  was  always  small  enough  to  ensure  sonic  flow.  The  orifice  calibra¬ 
tion  with  C2H2  was  used  for  the  C2D2  by  making  a  small  correction  for  the  mass 
difference.  Comparison  of  mass  spectra  from  C2H2  and  C2D2  flames  showed  a 
general  one-to-one  correspondence  for  the  peaks  in  the  same  general  pattern  of 
relative  intensities;  the  deuterated  ion  peaks  were  shifted  to  higher  mass  by 
a  difference  that  generally  increased  for  larger  ions. 

The  initial  mass  assignments  were  based  on  the  distance  in  the  recorded 
spectra  between  the  wel 1-characterized  C3H3+  (and  K+:  39  amu)  and  C13H*+  (165 

amu)  ion  peaks.  Although  this  did  not  give  the  correct  masses  elsewhere,  it 
provided  a  common  mass  scale  for  which  a  mass  correction  function  could  be 
developed.  It  also  provided  an  approximation  to  the  true  mass  scale  for  com¬ 
paring  the  isotopic  peaks.  For  example,  any  peaks  that  coincide  on  this  scale 
must  be  identical  in  mass.  The  strategy  was  to  correct  the  mass  scale  for 
known  masses,  e. g. ,  alkali  metals,  apply  the  correcting  equation  to  the  "raw" 
spectra,  and  then  compare  the  protonated  and  deuterated  flame  spectra.  The 
hydrogen  content  of  an  ion  is  given  by  the  shift  upon  deuteration,  while  the 
remaining  mass  indicates  the  number  of  carbon  atoms  and  the  possible  presence 
of  oxygen,  the  only  other  element  in  the  flame. 

Mass  assignments  were  made  to  be  consistent  with:  low  mass  hydrocarbon 
ions  identified  in  previous  work, 15  alkali  metal  and  lead  ions;  and  several 
integer  mass  relationships,  #1,  between  the  ions  observed  in  the  two  isotopic 
flames.  For  adjacent  species,  C,0»Hi  +  and  C, '0y'Hi/  +  (and  the  deuterated 
analogs) : 


AM*  *  12(X'-X)  +  16 (Y'-Y)  ♦  (Z'-Z) 


protonated 
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AM,  =  12<X'-X)  +  16(Y'-Y>  +  2(Z'-Z) 

AM„.  =  Z 
1h'.'  =  Z' 


deuterated 
inter-isotope 
inter- isotope 


The  analysis  was  complicated  by  the  low  mass  spectrometer  resolution  used 
for  high  throughput:  single  ion  peaks  had  widths  on  the  order  of  2  amu,  in¬ 
creasing  to  4  amu  at  higher  masses.  Mass  peaks  observed  at  every  carbon  num¬ 
ber  from  3  to  45  (the  high  mass  limit  of  the  instrument),  appeared  as  single 
peaks  in  the  protonated  spectra,  but  many  became  multiple  peaks  upon  deutera- 
tion. 


A  change  in  the  number  of  hydrogen  atoms  in  adjacent  peaks  is  indicated 
by  a  consistent  set  of  mass  changes.  For  example,  if  adjacent  peaks  differed 
by  CH2  then  AM.  =  14,  AM.  =  16,  and  AM..7  -  AMH»  *  2.  The  restrictions  on  the 
values  of  these  separations  provide  a  method  to  monitor  the  changing  hydrogen 
content  of  the  ions  throughout  the  mass  spectrum,  one  carbon  number  at  a  time. 
Ambiguities  were  resolved  through  considerat ion  of  the  effect  of  multiple, 
overlapping  mass  peaks.  An  example  is  provided  by  the  peaks  corresponding  to 
C35  and  C36  which  were  simulated  as  sums  of  simple  Gaussians  in  Fig.  1.  Each 
of  the  ion  species  will  exhibit  a  significant  i3C  ion  peak  with  intensity, 
relative  to  the  pure  i2C  ion,  proportional  to  the  number  of  carbon  atoms, 
e.  g. ,  1.1%  per  carbon  atom.  Each  observed  peak  was  then  treated  as  the  sum  of 
contributions  from  ion  species  differing  in  the  number  of  hydrogen  (deuterium) 
atoms  and  their  respective  *3C  isotope  peaks.  The  symbols  represent  the  ex¬ 
perimental  profiles  of  the  adjacent  C35  and  C34  peaks  from  the  C2H2  flame 
(Fig.  la)  and  the  C2D2  flame  (Fig.  lb).  The  apparent  single  peaks  in  Fig.  la 
are  compared  with  the  obviously  multiple  peaks  in  the  C2D2  flame  (Fig.  lb). 

The  sparse  lines  represent  calculated  profiles  for  individual  ions  which  sum 
up  to  yield  the  heavy  line  matching  the  experimental  data.  The  following 
constraints  were  observed:  (1)  Individual  ion  peaks  (e.g.,  C35HV)  were  approx¬ 
imated  as  Gaussians  with  widths  (FWHM)  of  4  amu  and  adjustable  amplitudes,  I0. 
As  drawn  in  the  figure,  each  peak  also  includes  the  contribution  from  the  *3C 
isotopic  ions  at  1  amu  higher  mass  and  (from  the  number  of  carbon  atoms)  a 
statistically  calculated  intensity  of  0.4  I0.  (2)  The  mass  scale  (amu/chart 

distance)  had  the  same  value  for  both  Fig.  la  and  lb.  (3)  The  chemical  formu¬ 
las  of  all  the  individual  ions  were  identical  for  the  two  flames.  In  this 
case,  the  formulas  in  Fig.  la,  C3jH,,  CjjHt.  1,  C3sHy»2,  ,  C3^H,.|  and 

C3jH».2,  become  C3$Dy,  C33Dy.it  C3sDy.2,  C34DY,  C3gDv.i  and  C3fDy»2  in  Fig. 
lb.  The  relative  intensities  of  the  individual  ion  peaks  were  varied  to  obtain 
a  good  match  between  the  calculated  total  profiles  (heavy  lines)  and  the  ex¬ 
perimental  profiles  (symbols).  The  deduced  intensity  patterns  are  not  identi¬ 
cal  for  the  two  flames  (perhaps  due  to  somewhat  different  flame  conditions) 
yet  the  spectra  both  are  dominated  by  C3SHy,  C35HT.2,  C36Hy,  C36HV.2,  and  the 
corresponding  deuterated  ions,  with  intensities  of  y  greater  than  y+2  for  C35 
and  y  less  than  y+2  for  C34.  Small  contribut ions  from  ions  with  Hy.i  may  be 
indicated,  but  we  feel  unjustified  on  the  basis  of  the  poor  spectral  resolu¬ 
tion  to  accept  their  presence  as  confirmed.  For  all  carbon  numbers  above  Ci3, 
the  deuterated  peaks  were  considerably  broader  than  the  protonated  peaks. 
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This  analysis  indicates  that  the  peaks  include  at  least  two  different  ionic 
species  which  differ  by  two  hydrogen  atoms,  corresponding  to  2  amu  and  4  amu 
in  the  C2H2  and  C2D2  flames,  respectively.  Oxygenated  species  cannot  be  ruled 
out  as  also  being  present,  for  if  C1.jHv.2O  (with  the  same  mass  as  CiHv.2)  is 
deuterated,  the  new  ion  will  have  the  same  mass  as  C« DY .  This  would  only  show 
up  as  a  shift  of  relative  intensity  to  this  first  peak  coincident  with  C,0T. 
Observation  of  this  was  difficult  since  the  intensity  pattern  in  the  protonat- 
ed  peaks  could  not  be  discerned  directly. 

The  simulations  presented  in  Fig.  1  also  show  that  the  C35  and  C34  pairs 
differ  only  by  one  carbon  atom  and  have  the  same  hydrogen  contents.  Because 
of  the  switching  of  the  dominant  ion  (C3SHY  and  C3«Hv*2>  the  apparent  C3S  to 
C34  mass  difference  for  the  overall  peaks  appear  to  be  larger  than  12  amu.  If 
the  alternate  peaks  dominated  (e. g. ,  C36Hv.2  and  C37HY  dominating)  a  value 
less  than  12  amu  would  be  apparent.  This  "switching"  phenomenon  was  observed 
in  the  protonated  spectrum  and  further  supports  the  contribution  of  two  ion 
masses  to  the  apparent  ion  peaks,  and  it  provides  a  rationale  by  which  to 
judge  which  ion  peaks  exhibit  a  change  in  hydrogen  content.  Using  an  "aufbau" 
approach  we  were  able  to  follow  the  chemical  composition  change  from  the  well 
known  C13H,+  ion  to  the  high  mass  limit  of  the  spectrometer.  From  the  compo¬ 
sition  of  the  ions,  the  masses  were  determined  and  were  found  to  be  as  much  as 
5-10  amu  lower  than  would  be  determined  using  a  simple  linear  extrapolat ion  of 
the  low  mass  data  to  the  2  500  amu  range.  The  polynomials  used  to  correct  the 
initially  deduced  masses,  M0,  were,  for  the  mass  range  to  165  amu: 


M, Rut  ■  -2. 9983  +  1.1149M,  -  9.2304  x  10‘*  M„ 2  +  2.0143  x  10-*M, 3 

and  for  the  mass  range  above  165  amu: 

Mt.m.  *  6.1210  +  0. 9594M*  +  6.0489  x  10-5M«2  -  9.5347  x  10-«M0 3 

Our  best  estimates  for'  the  large  ion  identities  are  given  in  Table  II.  Sug¬ 
gested  structures  for  all  of  these  ions  are  presented  in  Table  III.  Alternate 
structures  are  possible  for  some  of  the  ions  but  given  the  C/H  ratio  and  the 
rules  for  organic  chemistry  there  is  not  as  much  freedom  in  choosing  a  struc¬ 
ture  as  one  might  initially  suppose.  The  largest  ion  observed,  C4SH17,  mass 
559,  is  still  small  compared  to  the  smallest  soot  particle  observed,  The  mo¬ 
lecular  ions  in  Table  III  would  most  probably  be  flat  plates  with  diameters  up 
to  about  1.0  nm  for  mass  559  compared  with  the  1.5  nm  diameter  of  the  smallest 
observed  soot  particle. 3 


C.  MASS  SPECTRA 


We  have  obtained  ion  spectra  for  acetylene-oxygen  flames  at  2.7  kPa,  50 
cm/s  unburned  gas  flow  for  0«  1.0,  1.5,  2.0,  2.25,  2.5,  2.6,  2.7,  2.8,  2.9, 
and  3.0.  This  spans  the  range  from  nonsooting  to  well  into  the  sooting  range; 
the  soot  threshold  for  this  flame  is  at  0  =  2.4.  Some  of  the  results  are 
presented  in  Fig.  2  for  a  flame  at  0=  2.5.  The  other  results  will  not  be 
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reported  until  we  have  calibrated  the  system  so  that  the  mass  spectrometer  ion 
current  can  be  corrected  to  ion  concentrat ion.  This  should  be  done  in  the 
next  reporting  period.  Interpretat ion  of  the  data  will  also  wait  for  comple¬ 
tion  of  the  calibration. 

In  the  above  spectra  all  of  the  ions  contain  an  odd  number  of  hydrogens. 
This  is  contrasted  with  the  predictions  of  Stein1®  for  this  flame.  Stein 
argues  from  thermodynamics  and  kinetics  that  the  equilibrium: 

BPAH-H+  *  BPAH-+  +  H*,  (1) 

where  BPPH-H+  is  a  protonated  benzenoid  polyaromatic  hydrocarbon  (many  of  the 
species  in  the  above  sequence)  and  BPAH*+  is  an  ionized  benzenoid  polyaromat¬ 
ic,  is  always  to  the  right-hand  side  in  this  flame  in  the  temperature  range 
above  1500-1800  K.  He  accounts  for  our  observations  by  arguing  that  the  ions 
in  a  volume  of  the  flame  cooled  by  the  sampling  probe  are  cooled  below  *'1500  K 
so  that  BPAH-+  is  irreversibly  converted  to  BPAH-H+  in  several  microseconds. 
Assuming  a  rate  constant  of  10*10  cm®  molecule*1  s*1  for  the  protonation  reac¬ 
tion  he  deduces  a  half-life  of  3  x  10*®  s.  While  we  do  not  know  the  rate  of 
cooling  in  our  particular  nozzles,  Hayhurst  and  Telfrrd19  in  an  extensive 
analysis  for  similar  nozzles  estimate  a  flow  time  through  the  nozzle  of  about 
5  x  10*7  s.  It  thus  seems  reasonable  to  conclude  that  the  ions  we  observe  are 
the  actual  ions  in  the  flame  and  that  the  equilibrium  discussed  by  Stein  does 
not  play  a  major  role  in  shuffling  the  ions.  Further  note  that  the  observed 
ions  with  odd  numbers  of  hydrogen  atoms  are  just  what  is  predicted  assuming 
C3H3+  as  the  initial  ion. 


D.  REACTION  HEATING  OF  SOOT  PARTICLES 

In  the  ionic  mechanism  of  soot  formation  it  is  assumed  that  the  large 
ions  and  very  small  charged  soot  particles  are  formed  by  the  growth  of  smaller 
ions  by  the  addition  of  small  neutral  species,  such  as  acetylene.  Others 
attempt  to  explain  the  presence  of  these  ions  and  small  charged  particles  by 
other  mechanisms.  Homann  has  proposed11  that  as  the  incipient  soot  particles 
grow,  the  exothermic  growth  reactions  deposit  energy  in  the  particles  more 
rapidly  than  they  can  dissipate  the  energy.  The  particle  temperature  thus 
exceeds  the  gas  temperature  so  the  particles  are  thermally  ionized.  Homann 
has  not  developed  a  quantitative  description  of  the  phenomena®®  ;  it  is  incum¬ 
bent  upon  us  to  analyze  this  alternate  mechanism. 

We  thus  have  made  quantitative  estimates  of  the  amount  of  reactive  heat¬ 
ing  expected  for  growing  species  in  a  0  =  3.0,  2.67  kPa  C2H2/02  premixed  flame 
(50  cm  s*1  unburned  gas  velocity).  The  model  selected  is  similar  to  that  used 
by  Millikan21  to  estimate  particle  temperatures.  We  first  describe  the  model 
and  quantitative  estimates  of  the  various  input  parameters,  and  then  discuss 
the  implications  on  particle  heating. 

Growth  of  species  from  molecular  size  (150  amu)  to  particle  size  (10s 
amu)  is  considered.  P„,  designates  a  particle  which  has  undergone  n  reactive 
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collisions.  The  particle  is  heated  in  the  growth  reaction  with  species  G, 


P 


n-1 


<T0> 


k  <n) 

G  - > 


where  kft(n)  is  the  reaction  rate  constant  for  the  process  creating  P„.  The 
initial  temperature  of  the  precursor,  Ph.t,  is  T0,  and  the  entire  exother- 
micity,  q$,  of  the  reaction  is  assumed  to  heat  particle  Ph  of  heat  capacity, 
Cn* 


Particle  cooling  occurs  either  radiatively  or  by  collisions  with  the 
flame  species,  M,  as  described  by  the  process: 


k  (n) 

P  <T.)  +  M  — - - > 

n  1 


+  M 


where  qM  is  the  energy  transferred  to  the  bath  gas,  and  kM(n)  represents  the 
collision  rate  constant  of  M  with  P„. 

Several  properties  of  the  low  pressure  C2H2/02  flames  lead  to  simplifi¬ 
cation  of  the  model.  Relative  to  the  range  of  particle  sizes  considered  here 
(10  nm  diam),  the  gas  mean  free  path  is  very  large  (typically  on  the  order  of 
10  pm)  compared  to  the  particle  diameter  so  molecular  heat  transfer  processes 
dominate.  For  simplicity  we  neglect  radiative  cooling  of  the  hot  particles  so 
that  an  upper  limit  to  particle  heating  is  obtained. 

The  reactive  heating  rate  is  approximated  as  the  enthalpy  of  the  growth 
reaction,  AH*,  times  the  rate  of  reactive  collisions  of  the  growth  species  G 
with  the  particle:  AH$  k&CGl.  The  rate  constant,  k&,  can  be  written  as  the 
product  f&Z&CG]  where  Z*  is  the  collision  rate  of  G  with  the  particle  and  f& 
is  the  fraction  of  collisions  which  are  reactive.  Similarly,  the  rate  of  heat 
transfer  to  the  flame  gases  is  approximated  as  qMZMCM]  where  qM  is  the  heat 
transfer  per  collision  between  the  particle  and  M,  and  ZM  is  the  collision 
rr. te.  qM  is  given  by: 


qx  *  °  Cm  AT 


(1) 


where  a  is  the  accommodation  coefficient  for  energy  transfer  from  the  particle 
to  M,  CM  is  the  heat  capacity  of  ft,  and  AT  represents  temperature  difference 
between  the  particle  and  the  flame,  T*  -  T» . 


The  competition  of  these  two  processes  will  result  in  particle  heating  to 
a  steady-state  temperature  determined  by  the  condition: 

AH$f gZg  CG]  +  qMZM  CM]  *  0. 
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The  particle  temperature  is  then  given  by: 


T» («)  represents  the  "steady-state”  temperature  of  the  particles;  the  particle 
"temperature"  will  initially  be  a  function  of  time. 

We  first  consider  the  steady-state  particle  temperatures  for  various 
"particle"  sizes  from  molecular  size  (ISO  amu)  to  10s  amu,  corresponding  to  a 
spherical  particle  diameter,  d» ,  of  6  nm  (density  =1.5).  It  will  be  shown 
that  to  a  first  approximation,  T» (•)  is  independent  of  the  particle  size. 

Thus  the  particle  temperature  can  be  considered  independent  of  its  growth 
history  once  “steady-state"  is  reached.  The  denominator  of  the  right  hand  term 
of  Eq.  (2)  is: 

ac  Z  [Ml  =  y  B.C.Z.X.tMl  (3) 

MM  L  llii 

i 

where  the  sum  is  taken  over  the  major  flame  components  with  mole  fractions  X, , 
and  CM]  is  the  total  gas  concentration.  For  the  0  =  3. 0  flame,  the  major 
species  are  CO,  H2  and  C2H2  in  approximate  ratio5*9  of  6:3:1  and  CM]  =  1  x 
1017  molecules  cm*3  at  2000  K  and  2.7  kPa.  The  collision  number  Z,  from  ki¬ 
netic  theory  is: 


where  p,  is  the  reduced  mass  of  the  particle/reactant  pair.  For  large  parti¬ 
cles  this  smoothly  goes  over  into  the  rate  of  molecules  i  colliding  with  the 
surface  of  a  spherical  particle, 


Z  . 

l 


(5) 


Even  for  a  particle  of  mass  150,  use  of  mt  rather  than  p,  results  in  less  than 
10*  error  in  the  collision  number,  and  the  variation  in  di  for  H2,  CO,  and 
C2H2  (0.2-0.  A  nm)32  is  small  relative  to  d»  ('•  1  nm).  Assuming  a  to  be  inde¬ 
pendent  of  the  gas,  we  calculate  an  effective  value  for  Eq.  (3)  at  2000  K  in 
terms  of  the  CO  diameter,  dc •  which  will  be  essentially  valid  for  all  particle 
sizes.  The  heat  capacities  used  were  estimated  to  be  (kJ  K*1  mole*1)  26,  29, 
and  75  for  H2,  CO,  and  C2H2,  respectively.23  Likewise  the  reaction  collision 
number  Z»  can  be  simplified.  Equation  (2)  is  now  (employing  m&  '  26  amu  for 
acetylene,  the  only  real  growth  candidate) 2**2‘ 
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TP*  tf 


19  fGW  V  CM3 

""  a  <dp  +  dCQ>2  CM3 


(6) 


where  is  in  kJ  mole-i  and  the  diameters  are  in  angstroms.  Assuming  d&  * 
dco  <  d» ,  or  when  d»  is  very  large,  the  size  dependence  (except  any  hidden  in 
a)  drops  out  and 


19<AH  'k j 
G 


mole-1) 


Vg 


(7) 


This  equation  holds  equally  well  for  molecules  and  particles  from  150-10*  amu 
in  a  0  =  3.0  flame.  Now,  in  order  to  evaluate  T, ,  we  need  to  estimate  the 
heat  of  reaction  which  we  have  assumed  involves  the  major  hydrocarbon  species, 
C2H2.  We  note  that  mass  spectrometric  measurements  in  low  pressure  C2H2/02 
flames  have  shown  that  the  next  most  abundant  hydrocarbon  species,  C4H2,  is  an 
order  of  magnitude  lower  in  concentrat ion. * » 9 

Mature  soot  particles  show  considerable  graphitic  nature  while  small, 
young  particles  have  higher  H/C  ratios  which  decrease  as  the  particle  grows. 

We  estimate  the  exothermicity  of  the  growth  reactions  by  two  overall  process¬ 
es. 


C  2H  2 

1/4C6H6  +  C2H2 


-♦  2C  (s)  +  H2 

-*•  l/4CnH10  +  1  /2HZ 


AH2?b  =  -227  kJ  mole- 1 
AH2,8  =  "178  kJ  mole*1 


Using  available  heat  capacity  data, 23,27  at  2000  K,  the  processes  have  exo- 
thermicities  of  about  210  and  165  kJ  mole-i,  respectively.  Thus  Eq.  (7)  be¬ 
comes: 

400f 

T_  S  T_  +  - —  (8) 

P  F  o 

for  the  more  exothermic  process. 

The  coefficients,  f$  and  a,  remain  the  major  uncertaint ies  as  Millikan 
found. 2 l 

Energy  transfer  in  reacting  molecular  systems  has  been  investigated  both 
for  thermal  and  chemically  activated  systems.  The  efficiency  of  energy  trans¬ 
fer  is  expressed  in  terms  of  an  average  energy  step  based  on  a  model  of  energy 
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transfer  probability,  A  number  of  investigations  have  shown  that  the  mean 
energy  transferred  per  collision  decreases  with  increasing  temperature. 28-30 
Particularly  apropos  to  the  accommodation  coefficient  a  for  present  considera¬ 
tions  are  experiments  by  Barker  and  coworkers2**31  on  vibrational  energy 
transfer  from  highly  excited  ground  electronic  state  azulene  produced  in  in¬ 
ternal  conversion  of  the  laser  excited  molecule.  This  mimics  the  exothermic 
growth  process  in  that  the  vibrational  temperature  of  the  large  molecule  is 
much  higher  than  that  of  the  bath  gas.  Using  different  wavelengths  for  ex¬ 
citation,  the  dependence  of  the  transferred  energy  on  the  azulene  energy  con¬ 
tent  was  probed.31  The  average  collisional  energy  transfer  to  CO  was  about 
375  cm*1  for  azulene  excitation  of  30,600  cm*1,  and  it  was  150  cm*1  for  17,500 
cm*1  excitation.  These  values  correspond  to  a  »  0.2  and  0.1  for  CO  in  the 
high  and  low  energy  excitations,  respectively.  The  corresponding  values 
(energy  transferred  and  a)  for  H2  (about  the  same  heat  capacity  as  CO)  were 
roughly  half  of  these.  Barker  and  Golden  2*  found  roughly  a  T*°*  5  dependence 
for  the  average  energy  transferred  to  the  diatomic  N2  in  the  azulene  system. 

In  the  range  300-525  K  it  was  about  295  cm*1  at  room  temperature,  and  it  drop¬ 
ped  to  100  cm*1  from  400  to  525  K  in  good  agreement  with  the  observations  of 
Brown  et  al. 30  for  N2  in  ethyl  acetate  from  300-800  K.  Also,  Heymann  et  al. 32 
found  the  average  energy  transferred  to  CQ  from  cycloheptatriene  to  exhibit  a 
roughly  T*°*  5  dependence  while  other  monatomic  and  diatomics  exhibited  some¬ 
what  lower  temperature  dependencies.  Considering  the  major  post  flame  compo¬ 
nent  CO,  we  correct  the  conservative  a  =  0. 1  to  flame  temperatures  by  scaling 
by  the  inverse  square  root  of  the  temperature  which  gives  a  *  0.04. 

Gas-surface  accommodation  coefficients,  fg.,  are  generally  fairly  high  at 
room  temperature  (0.3-1)  and  also  have  been  observed  to  decrease  with  increas¬ 
ing  temperatures.  In  a  study  of  gas-surface  single-collision  energy  transfer, 
Rabinovitch  and  coworkers23  found  a  vibrational  accommodation  coefficient  of 
0.2  for  cyclobutene  on  quartz  at  975  K.  We  assume  that  a  for  the  particles  is 
at  least  as  large  as  the  molecular  value  of  0.04  estimated  above. 

The  particle  growth  rate  is  obtained  from  the  electron  microscopic  stud¬ 
ies  by  Howard  and  coworkers. 4  From  their  results  they  calculated  surface 
growth  rates  for  soot  particles  in  the  "standard"  acetylene-oxygen  flame  (P  = 
2.7  kPa,  u  =  50  cm/s)  below  about  2.5  pm/s.  This  corresponds  with  an  estimated 
value  of  f&  on  the  order  of  0.01  for  C2H2  as  the  growth  species.  This  com¬ 
pares  with  a  maximum  value  of  fg  =  0.001  (for  C2H2)  from  the  atmospheric  pres¬ 
sure  flame  measurements  of  Harris  and  Weiner. 24-24 

Substituting  a  =  0.04  and  fg  =  0.01  in  Eq.  (8)  gives  a  temperature  for  the 
particle  100  K  greater  than  the  gas  temperature.  If  we  increase  the  measured 
gas  temperature  by  100  K  and  calculate  the  total  ion  concentrat ion  using 
Saha’s  equation  the  results  reported  in  Fig.  3.  are  obtained.  The  small 
increase  in  equilibrium  concentrat ion  of  charged  particles  still  gives  a  theo¬ 
retical  concentrat ion  below  the  experimental  value.  The  decrease  in  the  mea¬ 
sured  value  with  distance  is  consistent  with  the  picture  that  the  ions  are 
formed  in  the  flame  front  and  then  disappear  by  ion-electron  recombinat ion. 

At  further  distances  from  the  burner  the  measured  and  calculated  charged 
particle  concentrations  may  agree  because  when  the  particles  become  larger 
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thermal  equilibrium  ionization  prevails.  This  is  completely  different  from 
the  ionization  we  are  concerned  about  in  the  earlier  part  of  the  flame;  unfor¬ 
tunately  the  two  parts  of  the  flame  are  often  confused. 

The  results  reported  in  Fig.  3  do  not  support  Homann’ s  proposal. **  We 
are  not  as  confident  as  we  would  like  to  be  however  with  the  values  of  a  and 
f*  so  have  held  up  publication  until  we  can  find  more  reliable  values. 


III.  PUBLICATIONS 


The  following  papers  were  published  during  this  report  period: 

1.  "Ionization  and  Soot  Formation  in  Premixed  Flames,"  D.G.  Keil,  R. J. 

Gill,  D. B.  Olson,  and  H. F.  Calcote,  Twentieth  Symposium  (International) 
on  Combustion  (The  Combustion  Institute,  Pittsburgh,  1985)  p.  1129. 

2.  "The  Effect  of  Temperature  on  Soot  Formation  in  Premixed  Flames,"  D. B. 
Olson  and  S.  Madronich,  Combust.  Flame  60,  203  (1985). 

The  following  manuscript  is  in  preparation: 

1.  "Langmuir  and  Thermocouple  Probe  Measurements  in  Flames,"  D.G.  Keil, 

R. J.  Gill,  D. B.  Olson,  and  H.  F.  Calcote.. 


IV.  PERSONNEL 

In  addition  to  the  authors  the  following  personnel  made  significant  con¬ 
tributions  to  this  program: 

D.  B.  Olson,  Physical  Chemist 

H.  Rothschild,  Librarian,  Technical  Editor 
(Ms.  Rothschild  also  assisted  with  data  reduction.) 

E. L.  Stokes,  Technical  Typist 


V.  TECHNICAL  INTERACTIONS 


Technical  interactions  with  the  scientific  community  have  taken  several 
forms;  the  most  important  are  presentations  at  scientific  meetings,  seminars 
and  workshops.  Proposal  and  manuscript  review  has  also  played  an  important 
role  during  this  report  period.  Olson  and  Calcote  visited  Dr.  Charles  Martel 
at  Wright  Field  to  discuss  our  work,  especially  the  effect  of  fuel  molecular 
structure  on  jet  engine  performance.  The  following  people  visited  AeroChem 
during  this  contract  to  discuss  the  mechanism  of  soot  formation  at  a  workshop 
sponsored  by  Army  Research  Office: 
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T.  Brabbs  and  E.  Lezberg  (NASA/Lewis  Research  Center) 

J.  Eyler  (University  of  Florida 

W.  Flower  (Sandia  National  Laboratories) 

M.  Frenklach  and  T.  Lester  (Louisiana  State  University) 

I.  Glassman,  A.  Gomez,  and  G.  Sidebotham  (Princeton  University) 

S.  Harris  (General  Motors  Research  Laboratories) 

R.  D.  Kern  (University  of  New  Orleans) 

U.  G.  Mallard  and  R.  Santoro  (National  Bureau  of  Standards) 

D.  Mann  (U. S.  Army  Research  Office) 

J.  Tishkoff  (Air  Force  Office  of  Scientific  Research) 

The  following  present at  ions  were  made: 

1.  "Soot  Formation  in  Flames"  Presented  by  H.  F.  Calcote  to  six  separate 
groups  on  a  "People  To  People"  tour  through  China  in  which  Calcote  was 
the  Delegation  Leader,  9-30  October  1984. 

2.  "The  Effect  of  Temperature  on  Soot  Formation  in  Premixed  Flames"  D.  B. 

Olson,  Eastern  States  Section:  The  Combustion  Institute,  Fall  Technical 
Meeting,  Clearwater,  FL,  December  1984. 

3.  "The  Role  of  Ions  and  Charged  Particles  in  Soot  Formation"  H.  F.  Calcote, 
Seminar  Rutgers  University,  Engineering  and  Aerospace  Dept.  10  April  1985 

4.  "An  Ionic  Mechanism  of  Soot  Formation  in  Flames"  Calcote,  invited 

presentation  at  the  Particle  Emission  Technology  Meeting",  Naval  Post 
graduate  School,  Monterey,  CA,  16-18  April  1985 


VI.  INVENTIONS  AND  PATENT  DISCLOSURES 

There  are  no  inventions  or  patent  disclosures  to  report  during  this  re¬ 
porting  period. 
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TABLE  I 

WHY  I  BELIEVE  IN  THE  IONIC  MECHANISM  OF  SPOT  FORMATION  IN  FLAMES 

1.  Fuel  effects:  the  propensity  to  form  soot  correlates  with  the  tendency  of 
the  fuel  to  produce  ions. 

2.  Chemical  additive  effects:  the  explanation  of  some  additive  effects  is 
contingent  upon  an  ionic  mechanism  of  incipient  soot  formation  -  not  to 
be  confused  with  ionic  additive  effects  on  soot  agglomeration. 

3.  Soot  formation  is  a  nonequi 1 ibrium  process  at  the  critical  equivalence 
ratio  for  soot  formation;  it  thus  needs  a  nonequilibrium  process  to 
drive  it — chemi-ionization  is  a  nonequilibrium  process. 

4.  Ion-molecule  reactions  to  explain  the  phenomena  are  easy  to  prescribe,  no 
ad  hoc  assumptions  are  required — compare  the  problems  with  free  radical 
reactions.  For  example: 

a.  The  ions  required  in  the  mechanism  have  been  observed  in  flam¬ 
es. 

b.  Ion-molecule  reactions  are  fast. 

c.  Isomerizat ion  of  ions  is  very  fast;  explains  the  transition 
from  a  linear  structure  to  a  cyclic  structure  -  a  major  problem 
with  free  radical  mechanisms. 

5.  The  concentrat ion  and  rates  of  ion  formation  seem  to  be  adequate  to 
account  for  the  appearance  of  soot — there  may  be  some  questions  here 
because  of  the  limited  number  of  tests  and  the  quality  of  the  soot  con¬ 
centration  data. 

6.  Dramatic  changes  occur  in  the  ions  at  the  soot  threshold,  consistent  with 
what  would  be  expected  for  an  ionic  mechanism. 


7. 


Variation  in  ion  composition  and  location  of  ion  profiles  with  respect  to 
soot  formation  is  consistent  with  an  ionic  mechanism  of  formation  of 
incipient  soot. 
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TABLE  II 

ION  IDENTITIES  OF  MAJOR  HEAVY  IONS  OBSERVED  IN  RICH  C2H2/02  FLAMES 


C 1 #H 1 1 

it  a* 

C 1 »H  1 1 

•  t3+ 

CaoHn 

.  13+ 

CaiHi 1 

.  13+ 

C  a  *H 1 3 

,1S+ 

C23H13 

.t5+ 

c24h13 

.1S+ 

c2Sh13 

.IS4" 

Ca*Hia 

.IS+ 

C37Hi3 

.13+ 

C2tH13 

.IS4" 

C  a  »W  1 3 

,1S+ 

C30H1S 

,17+ 

c31h15 

,17+ 

Csa^i j,j7+ 
G33H15, 174- 
C34H1S» 17+ 
CssHlS. 17+ 
C34H15» 17+ 
C37HlS.17+ 
C3»Ht5, 17+ 

Cat^is, t7+ 
C 40^17. l»+ 

C4iH17.lt"*" 
C  4  2^  i  7  »  1  t* 
C  4  3H  i  7  •  1 t+ 
C44H17  » lt+ 
C45Hi7.lt"*" 
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TABLE  III 


SUGGESTED  STRUCTURES  OF  OBSERVED  FLAME  IONS 


C18Hii+<227) 


CiaH13+  (229) 


ramiaC 

Ci9Hh+(239) 


jqmC, 

Cl  9Hi3f  (241) 


DIO 


DIO 


C2 O^i  (  2  5  3  ) 


DJl+J 


OIQ 


02(^13^  (255) 


opjo 

I0IDI 


C21Hn+  (265) 


C2iHi3+(267) 


C34H15*  (423) 


C34Hi7+  (425) 


xngr 

X^mqC 

xSsoy 

'ijonf 


C35H15  (435) 


C3  5H 1  7+  (437) 


C37H1  5+  (459) 


C37Hi7+(46  1) 


oro 


DIO 


OlOlQ 


DID 


DID 


OJDJPI  ^ 


dtoioiq 


DT+uD. 


DIOIOTOTO 


OIQ 


Dior  TO 


aormpin 


oiqlojo 


AMPLITUDE,  arbitrary  units  AMPLITUDE,  arbitrary  units 
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MASS  SHIFT  FROM  CggHy,  amu 


(a) 


MASS  SHIFT  FROM  amu 

(b) 

FIGURE  1  COMPARISONS  OF  CALCULATED  AND  EXPERIMENTAL 
ION  PEAKS.  C2H2  AND  C2D2  FLAMES 

Sparse  lines  -  simulated  peaks  for  species  indicated,  (relative 
intensities  in  parentheses) .  Heavy  lines  -  sum  of  calculated 
individual  peaks.  □  -  experimental  profiles  (see  text). 


